Spin models for orientational ordering of colloidal molecular crystals 
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Two-dimensional colloidal suspensions exposed to periodic external fields exhibit a variety of 
molecular crystalline phases. There two or more colloids assemble at lattice sites of potential minima 
to build new structural entities, referred to as molecules. Using the strength of the potential and 
the filling fraction as control parameter, phase transition to unconventional orientationally ordered 
states can be induced. We introduce an approach that focuses at the discrete set of orientational 
states relevant for the phase ordering. The orientationally ordered states are mapped to classical 
spin systems. We construct effective hamiltonians for dimeric and trimeric molecules on triangular 
lattices suitable for a statistical mechanics discussion. A mean-field analysis produces a rich phase 
behavior which is substantiated by Monte Carlo simulations. 



PACS numbers: 82.70.Dd, 64.70.Dv 

I. INTRODUCTION 

Soft materials comprised of colloidal particles undergo 
phase transitions from fluid to crystalline order just as 
their atomic counterparts 0, El 3- The intrinsic time 
and length scales of such suspensions offer the advan- 
tage to monitor the ordering phenomena by video mi- 
croscopy on a single particle level. Furthermore, the in- 
teractions can be tailored to a certain degree by con- 
trolling the surface charge of the colloids or changing 
the refractive index of the particles and the solvent or 
by screening electrostatic interaction by the addition of 
salt, etc. H, 0- Colloids thus constitute fascinating 
model systems to study material properties. For ex- 
ample, new crystalline structures without known atomic 
analogue have been found Q for mixtures of charged col- 
loids. 

In particular, confining geometries allow to study 
phase transitions in two-dimensional systems where 
strong thermal fluctuations often become the relevant 
mechanism for the melting of order. For example, close to 
melting a universal relation betwen the elastic constants 
is predicted in the framework of the Kosterlitz-Thouless- 
Halperin-Nelson- Young theory |1, Pp, Ho| as has been re- 
cently observed experimentally |ll| . Even richer phase 
behavior is expected if quasi two-dimensional colloidal 
systems are exposed to external potentials as can be re- 
alized by interfering laser beams 0, 0, periodic 
pinning arrays |15| . imprint or stamping techniques |16| . 
For example, in the case of one-dimensional troughs ex- 
otic phases such as the locked smectic phase or a floating 
solid have been predicted [l7L Il8| and partially observed 
in simulations 19, 20] and experiments [21|. 

For three or more laser beams periodic patterns with 
two modulation directions can be achieved |22j| . in par- 
ticular, square and triangular lattices have been stud- 
ied IH, 0, ■ The colloidal particles then accumu- 
late in the potential minima thus imposing a modulation 



of the colloid density. The competition between the at- 
traction to the potential minima and the mutual repul- 
sion of the particles can be conveniently controlled by 
adjusting the filling factor, i.e., the number of colloidal 
particles per potential minimum. For instance, by chang- 
ing the angles of incidence the lattice constant of the pe- 
riodic modulation can be varied continuously. Since the 
electrostatic interaction is usually screened by counteri- 
ons in the solvent, the lattice spacing also sensitively de- 
termines how strongly particles interact. A second route 
to affect the balance between external compression and 
internal repulsion is by changing the intensities of the 
laser beams. Third the salt concentration of the solvent 
is a suitable control parameter to tune the range of the 
repulsive interaction, i.e. the Debye screening length. 

This paper is motivated by recent experimental studies 
on two-dimensional colloidal system exposed to a poten- 
tial of a triangular lattice symmetry with a filling factor 
three . The density of the colloids has been selected 
in a regime below the spontaneous freezing transition. A 
continuous increase of the strength of the laser field has 
induced a sequence of qualitative changes of order. At 
low intensities of the lasers, corresponding to a potential 
depth of say Vq ~ fc^T, the colloids exhibit the typical 
response of a fluid to an external modulation, i.e., a small 
periodic component is superimposed to the homogeneous 
background density. Enhancing the strength of the laser 
potential by an order of magnitude the response becomes 
highly nonlinear; groups of three colloidal particles clus- 
ter in the close vicinity of each potential minimum. How- 
ever, a significant number of defects - groups of two or 
four colloidal particles - is still present in the system. 
Since the filling factor has been adjusted to three, the 
number of four-groups equals the one of two-groups. In 
this regime the barriers separating minima become high 
and correspondingly single particle diffusion is strongly 
suppressed H3 , similar to a solid. Although no thermo- 
dynamic phase transition has occurred the system is re- 
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ferred to as a molecular crystal. A further increase of the 
external potential results in a freeze-out of the defects, 
i.e., each lattice site is occupied by precisely three col- 
loids. The confinement is rather strong, balanced only by 
the mutual electrostatic repulsion, leading to triangular 
arrangements of the three-groups aligned with the lattice. 
Interestingly, the majority of three-groups is oriented in 
the same way with a small number of groups adopting 
a mirrored configuration, see Fig. ^ Using a suitable 
order parameter, orientational correlations demonstrate 
the long-range orientational order. At still higher laser 
potentials the orientation of the triangles is uncorrelated 
at large distances, indicating a phase transition to an 
orientationally disordered phase. In Ref. [22j this subse- 
quent melting has been referred to as reentrant transi- 
tion, see Fig. ^ 
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FIG. 1: Contour plots of the lateral density distribution (left) 
and the averaged local particle density (right) for different 
light potentials: V = 60 k B T (top) and V = 110 k B T (bot- 
tom). In the first case the majority of three-groups of colloids 
is oriented in the same way, whereas in the later case the ori- 
entation of the triangles is uncorrelated at large distances. At 
some intermediate Vb a phase transition occurs. The horizon- 
tal/vertical axes are x and y, respectively. All units are in 
fim. Figure is reproduced from a part of Figure 2 of Brunner 
and Bechinger Ref. 122^1 . 

The possibility to induce phase transitions by peri- 
odic laser patterns has stimulated theoretical studies in 
two-dimensional colloidal systems. Reichhardt and Ol- 
son have performed extensive Langevin simulations and 
have identified ground states for a series of filling factors 
on square and triangular lattices |23| . In particular, they 
have also observed the orientational ordered state cor- 
responding to the experiment discussed above. For non- 
integer filling factors the molecular crystals are composed 



of supercells, i.e., several types of groups of colloids are 
distributed periodically on the optical lattice [2(j- By 
lumping particles to rigid molecules, explicit minimiza- 
tions of intermolecular interactions in the presence of the 
optical potential have revealed several ordered ground 
states [24j . Furthermore, an energy functional for dimeric 
molecules on a square lattice has been introduced and an 
analogy to Ising behavior has been identified [24[ . 

The goal of the present paper is to elaborate a theoreti- 
cal framework for colloidal molecular crystals in external 
fields built from a statistical mechanics perspective, as 
outlined in |26j. First, to deal with the complexity of 
the system we focus on the essential low-energy degrees 
of freedom of colloidal suspension. From the experimen- 
tal observation one infers that the groups of colloids are 
to be considered as rather rigid entities as soon as the 
laser potential exceeds several tens of UbT. The entities 
are referred to as molecules or dimers, trimers, etc. The 
key idea is to keep only discrete orientational states of 
the molecule as fluctuating variables, i.e., relevant for a 
statistical mechanics model. The separation of energy 
scales, e.g. the condensation or binding energy to build 
such molecules is much larger than the thermal scale, 
suggests to ignore processes connected with the breaking 
apart of the structure. The theoretical description we 
wish to develop is of phenomenological nature, since it 
does not explain the formation of molecules. In particu- 
lar, the regime of low laser intensities where the molecu- 
lar crystal exhibits a significant fraction of defects is not 
within the scope of our approach. The benefit of focusing 
on effective low-energy excitations is that one can derive 
simple models that are amenable to powerful methods of 
the statistical mechanics toolbox. In particular, it allows 
to identify the appropriate broken symmetry phases and 
to map the problem to related magnetic transitions for 
which a great wealth of knowledge has been developed. 

Specifically, in Section II we develop a theory for the 
experiment performed by Brunner and Bechinger |22| . 
We identify the relevant excitations of the system at the 
thermal scale as single trimer flips. A hamiltonian is 
constructed that couples nearest neighboring trimers via 
their respective orientations. The microscopic interac- 
tions, that are for example also responsible for the for- 
mation of the molecules, enter the theory only implicitly 
through three energy scales. The explicit link to the ex- 
perimental parameters, e.g., salt concentration, effective 
surface charge, etc., is deferred to the Appendix. It turns 
out that the statistical mechanics of trimers on a trian- 
gular lattice maps to an Ising model once the spectrum 
of single trimer excitations are evaluated. A comparison 
to experimental results corroborates our approach. 

In Section III we propose an experimental setup using 
a triangular lattice with a filling factor two. Then dimers 
are the composite objects and three orientational states 
per lattice site emerge. The corresponding hamiltonian 
can be understood as a natural extension of a three-states 
Potts model with the peculiarity that orientational pair 
interactions depend on the lattice orientation. We find 
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a rich scenario of four broken symmetry phases occupy- 
ing different regions in the phase diagram. In addition 
to "ferro-" and "antiferromagnetic" structures a herring- 
bone structure appears and a novel phase referred to as 
Japanese 6 in 1 (J6/1) is discovered. An analytical ap- 
proach based on a variational mean-field calculation is 
presented to explore the phase boundaries. Extensive 
Monte Carlo simulations are employed to corroborate our 
findings and we present a careful discussion of the differ- 
ent phases. 

In the Conclusion, we summarize our main results and 
put them in a broader perspective of experimental and 
theoretical solid state physics. 




(b) 



E 3 
E 2 
E^ 



S=1 



FIG. 2: The model system, (a) The triangular lattice of the 
external field and the two discrete orientational states of the 
trimers. (b) Three interaction energies for the orientational 
configurations of neighboring trimers. 



II. TRIMERS 



For densities of the colloidal particles adjusted such 
that there are three particles per potential minimum, 
trimers will form spontaneously at sufficiently strong 
laser intensity. Then the low-energy effective degrees of 
freedom are given by the discrete orientational states of 
the composite object. Due to the symmetry of the optical 
lattice the trimers constitute equilateral triangles which 
are aligned with the triangular lattice. From Fig. [21 one 
infers two energetically equivalent orientational states of 
the trimer corresponding to the optimal balance between 
the compression of the constituent colloids due to the 
laser potential and their mutual screened Coulomb repul- 
sion. The residual interaction between trimers depends 
on their relative orientation as well as on their relative 
lattice position. Since the residual interaction is solely 
due to a screened Coulomb interaction the magnitude 
of this energy decreases rapidly with increasing distance. 
For typical experimental setups the screening length is 
an order of magnitude smaller than the lattice constant 
[22|. which allows to restrict the residual interaction to 
nearest neighbors. For each neighboring pair of trimers 
there are four possible geometrical configurations, two of 
which are mirror images of each other, see Fig. [21 These 
configurations introduce three energy scales, with generic 
ordering E\ < E2 < E3 . These scales are functions of the 
screening length, the strength of the laser potential, as 
well as the effective charges of the constituent colloids. 
Appendix lAl provides a derivation of the explicit relation 
and supplements a careful discussion of the concept of 
composite objects. 

The considerations of the previous paragraphs allow 
to define the model in terms of a statistical mechan- 
ics problem. Each trimer i is attributed a "spin state" 
Si = ±1, the total number of orientational configurations 
being given by 2 N for a triangular lattice of N sites. The 
phase behavior of the trimer system is encoded in a "spin 
hamiltonian" , which reduces to a sum over local interac- 
tion energies. The trimer system exhibits a peculiarity in 
the sense that these energies are dependent on the bond 
vector, i.e., the direction of the connecting lattice vector, 
in addition to the two spin orientations of the neighbor- 
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FIG. 3: The two classes of equivalent bond vectors in the 
trimer system. Note that the classes have the same symmetry 
as the corresponding trimer state, i.e., for Si = a. 



ing spins. Explicitly, a simultaneous flip of both spin is 
not a symmetry operation of the hamiltonian, however, 
an additional rotation of the corresponding bond vector 
by 60° restores the original configuration. Consequently, 
a rotation of the bond vector by 120° with identical spin 
states results in the same interaction energy. Therefore, 
similar to the two orientational states of a trimer, there 
are only two inequivalent bond vectors, see Fig. [31 The 
hamiltonian can be expressed in the following form, 



a (ij)a 



(1) 



where a = ±1 denotes the orientation of the bond vec- 
tor, and the symbol (ij) a indicates the neighboring pair 
of spins with directed bond vector a. Hence, each bond 
appears precisely once in the entire sum. Reversing the 
order of the lattice sites implies an inversion of the bond 
vector (ij)a = (ji)—a- Explicitly, the local energy func- 
tional read 
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(2) 



Note that due to the above explained dependence on the 
orientation of the bond vector, the corresponding energy 
"matrices" are not symmetric. 

A global shift of the energy scale does not affect the 
phase diagram, leaving two independent energy differ- 
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ences. Introducing the thermal scale ksT allows to con- 
struct two dimensionless parameters, and naively one 
would expect a two-dimensional phase diagram. How- 
ever, we shall show now that due to geometrical con- 
straints there is a one-to-one correspondence of the ex- 
citation spectrum of the colloidal trimers and an Ising 
model on a triangular lattice. Since the interaction en- 
ergy is additive, it is sufficient to consider the case of 
three collinear trimers — i.e., considering the smallest 
neighborhood of trimers with representatives of both in- 
equivalent bond orientations. If we fix the central trimcr 
there are 2 2 possible configurations of the two neighbors, 
all of them are depicted in Fig. 01(a). Then, flipping the 
central trimer yields the energy changes ±AE, where 
we have defined AE = E\ + E3 — 2E2 . Although unex- 
pected due to the bond-oriented local interaction ener- 
gies, the trimer system reduces to an Ising model, i.e., 
with spins allowing for two distinct orientations and a 
single energy scale AE for excitations, see Fig. 0] (b). 
As a result of these considerations, the hamiltonian of 
Eq. 0} is equivalent to 
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H 



J ^2 SiSj 



Ho, 



(3) 



where ^2uj\ indicates the sum over all pairs of neigh- 
boring spins, irrespective of the orientation of the pair 
with respect to the lattice. The exchange energy J sets 
the scale for local excitations and is related to the three 
interaction energies by 

J=^AE=^(E 1 -2E 2 +E 3 ). (4) 

The additive constant can be evaluated to H.q = 
— jN(Ei + 2E 2 + £3). Although the energy scales pos- 
sess the generic ordering E-y < E 2 < E 3 , the sign of the 
exchange coupling J can attain in principle both nega- 
tive and positive values. For two-dimensional colloidal 
systems interacting via strongly screened Coulomb in- 
teraction explicit evaluation yields positive values, i.e., 
" ferromagnetic" coupling. For systems with different mi- 
croscopic interactions, e.g., with paramagnetic particles, 
an "anti-ferromagnetic" phase could possibly be realized. 

Let us emphasize that the mapping of the trimer hamil- 
tonian to an Ising model is different from the 1:1 corre- 
spondence of a lattice gas model for a binary alloy to 
an Ising system. There the nearest-neighbor interac- 
tion energy is parametrized also by three energy scales 
£aa j £ab , £bb corresponding to a pair of two neighbor- 
ing A atoms, a mixed pair of an A and a B atom, and 
two B atoms. For the mixed pair it is irrelevant whether 
atom A or B is left /right or up/down. For the trimer 
problem the configurations where trimers face each other 
or are back-to-back are energetically rather different, see 
Figure |3(b). 

For a two-dimensional Ising model an exact solution is 
available for various lattice symmetries |27| . On a trian- 
gular lattice the transition between the "ferromagnetic" 
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FIG. 4: (a) Schematic representation of all possible situations 
for a trimer flip with the corresponding energy changes; AE — 
Ei + Es — 2E2. (b) The same as in (a) but for spins of the 
spin-1/2 Ising hamiltonian Tt = — J^sUj) SiSj. 



and "paramagnetic" phase occurs at ksT c = 4J/ln3. In 
the experiments of Ref. , the temperature is kept fixed 
and the intensity of the external potential is varied. This 
variation is accompanied by the change of the linear ex- 
tensions of trimers, see Appendix^] which further affects 
the strength of the exchange energy J. Above some crit- 
ical laser intensity the exchange energy becomes smaller 
than the critical one, 



J < J c = k B T 



In 3 



(5) 



and the orientational order of the trimer system is lost. 

Employing experimental parameters from Ref. |22j . 
1/k = 570 nm for the inverse screening length, £q = 
3.4 x 10 4 ksT for the strength of the screened Coulomb 
interaction, and T — 295 K as room temperature (see also 
Appendix we have determined the phase diagram in 
the (Vo/ksT, (Ka)^ 1 ) plane and compared it to the ex- 
perimental results. Brunner and Bechinger report that 
for Vq = 60 fcsT the trimer system (realized for a lattice 
constant a = 11.5 /im) is orientationally ordered, whereas 
it is disordered for Vq = 110 ksT. The critical strength of 
the laser potential found in our calculations is in striking 
agreement with these observations, i.e., Vq c — 78.6 fcgT, 
see also Fig. A more systematic study, similar to what 
has been performed in the case of ID periodic laser po- 
tentials |28l |. would be highly desirable, in particular, in 
the regime of larger lattice constants and smaller screen- 
ing lengths where the shape of the phase boundary is 
more sensitive to control parameters. 

We have also compared our theory to Langevi n dy nam- 
ics simulation results of Rcichhardt and Olson |23j . The 
phase diagram is exhibited in Fig.|S]in the (Vo/£o, T/T® n ) 
plane. Although the simulation has been performed for 
dimers on a square lattice, the basic mechanism of the 
orientational melting - and thus the corresponding be- 
havior of the phase boundaries - is robust with respect 
to the symmetry of the composite objects and the un- 
derlying lattice. Note that the parameter Vq/Eq is by 
two orders of magnitude smaller than in the correspond- 
ing experimental setup; nevertheless our analytical solu- 
tion captures the main features of the phase boundary, 
i.e. shape and parameter range. The remaining phase 
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FIG. 5: Phase diagram of the colloidal trimer system as a 
function of the strength of the external potential (in units of 
IsbT where T = 295 K) and colloid-colloid interaction char- 
acterized by a product of the inverse Debye length and the 
lattice constant (solid line). The two crosses correspond to 
the experimental observations reported in Ref. |22T| . Insets: 
averaged local particle densities as represented in Figs. 2 (f) 
and (h) of Ref. fH . 

boundaries presented in the paper of Reichhardt and Ol- 
son are connected to the fission of the composite objects 
and are beyond the scope of our model. 



to the uniaxial deformation the exchange energy differs 
from the remaining two, J a and respectively. The 
system is still represented by a two-dimensional Ising 
hamiltonian, however, with anisotropic exchange ener- 
gies. This system again allows for an exact solution, in 
particular, the critical temperature is now determined by 
the equation [2?| 

e -4J a /k B T +2e -4J b /k B T = 1 _ (g) 

In conclusion, a small uniaxial strain leads only to a small 
shift of the critical temperature but does not lead to any 
qualitatively new phase behavior. 
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FIG. 7: Schematic representation of the lattice with one 
trimer and his neighbors for isotropic triangular lattice (mid- 
dle) and uniaxially deformed ones; below the angle between 
the horizontal bond and the non-horizontal one. Due to the 
symmetry, the exchange coupling in the horizontal direction 
differs from the one in the remaining directions. 
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FIG. 6: Phase diagram of the colloidal trimer system as a 
function of the reduced strength of the external potential and 
temperature. Our model describes part of the phase diagram 
corresponding to the orientational melting - solid line. Sym- 
bols correspond to the results from the simulation of Reich- 
hardt and Olson, see Ref. ^ Fig. 4.. 

Let us comment on the possibility to impose a small 
uniaxial strain on the triangular lattice by varying the an- 
gles of the incident laser beams. For small deformations 
the trimers will merely stretch leaving the system still 
with two orientational states per lattice site, see Fig. [7| 
Since the distances between pairs of trimers are now dif- 
ferent in different directions, the pair interactions become 
direction dependent. The analysis for the excitations en- 
ergies is now valid separately for each direction of the 
three co-linear trimers. For the direction perpendicular 



III. DIMERS 

In this Section we generalize the concepts introduced 
for trimers and discuss the possibility of dimeric molecu- 
lar crystals and their corresponding phase diagram. In- 
clining the incident laser beams with respect to the plane 
of the two-dimensional colloidal system one can adjust 
the lattice constant of the laser potential. Then, by 
appropriate matching of the lattice constant to the col- 
loidal particle density, one can achieve an average of two 
particles per potential minimum. For sufficiently large 
laser intensity, one expects that a defect-free structure 
is formed, i.e., the composite objects of the system are 
dimers. The symmetry of the lattice allows now for three 
different orientational states per lattice site denoted by 
o t = 1,2,3, see Fig. El (a). 



A. The dimer hamiltonian 

As in the case of trimers, the interaction of neighboring 
dimers is the origin of orientational ordering, whereas the 
balance of the internal repulsion and the laser pressure 
merely yields the binding energy of the dimers. From 
Fig. [S] (b) one infers that there are four interaction ener- 
gies for the 3 2 configurations of neighboring dimers with 
generic ordering Ei < E 2 < E 3 < E±. Explicit formula; 
for the relation to experimental parameters are presented 
in Appendix [S] The dependence of the interaction of 
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FIG. 8: The model system, (a) The triangular lattice of the 
external field and the three discrete orientational states of 
the dimers. (b) Four interaction energies for the orientational 
configurations of neighboring dimers. (c) The three classes 
of equivalent bond vectors in the dimer system. The groups 
have the same symmetry as the corresponding dimer state, 
i.e., for Gi = q. 



2 and from a = 2 to a = 3 dimer states have to be 
permuted cyclically, 1 — ► 2 — ► 3 — ► 1. 

As noted in [26| . the four interaction energies do not 
enter the phase diagram independently, but only two lin- 
ear combinations of these are relevant for the spectrum 
of excitations. Introducing the thermal scale fc^T, a two- 
dimensional phase diagram characterizes the complete 
orientational ordering scenario of the dimer system. Here 
we provide the full chain of arguments leading to this re- 
duction of control parameters. First it is convenient to 
sum over the three inequivalent lattice orientations sep- 
arately, i.e., to split the hamiltonian into TL = ^ Q Tt a 
with H a = J2(ij) c ,h( <J ii <J ji a )- Then the correspond- 
ing hamiltonian for each direction can be rewritten in 
a "spin nomenclature" as a generalized spin-l-Ising or 
generalized Potts model. Here we present the derivation 
for a = 1, the results for the remaining directions follow 
by cyclic permutations. The local energy functional can 
be represented as a sum of four simpler functionals 



h(<Ti,aj; 1) = y^^h r (ai,aj; 1), 



(9) 



dimers on the orientation of the bond vector follows from 
the symmetry of the composite objects, i.e., there are 
three inequivalent pairs of lattice orientations a — 1,2,3, 
depicted in Fig.|SI(c). Since each pair of equivalent lattice 
orientations is now collinear, reversing the order of sites 
in a directed bond (ij) a yields a directed bond within 
the same equivalence class. This implies that the local 
interaction energy is symmetric with respect to exchange 
of the dimers. 

The total hamiltonian of the two-dimensional dimer 
system corresponds to a sum over all interactions of 
neighboring dimer states. Denoting by (ij) a the pair 
of nearest neighbors i,j whose bond vector is parallel to 
the orientation of the dimer state a, the colloidal dimer 
hamiltonian can be expressed in the following form 



where the h r (ai, af, 1) = h r are given in terms of tables 



a=l (%j) a 



The direction-dependent local energy functionals 
h(<Ji,<jj;a) are given by: 
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Note that the h{<Ji, af, a) are symmetric with respect to 
the diagonal, and that they are related to each other 
through cyclic permutations, i.e., from a = 1 to a = 
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and each contribution is parameterized in terms of a sin- 
gle energy scale, 

K = -(Ei-E 2 ), 

L = E 3 -E 2l (11) 
M = -(E 4 -Ei -2E 3 + 2E 2 ). 

Due to the generic ordering of the interaction energies of 
colloidal dimers, the parameters K and L are expected 
to be positive, while M may attain both signs in gen- 
eral. The functionals allow for a direct interpretation: 
First, hi favors dimers to be in the same state, second, 
h 2 gives an extra contribution if a dimer is aligned with 
the bond vector connecting the dimers. Third, h 3 ac- 
counts for the high symmetry configuration, viz. both 
dimers aligned with the bond vector. Fourth, /14 is inde- 
pendent of the configuration and corresponds to global 
total shift of the energy. The functionals h r (oi , <jj ; a) are 
easily reexpressed in terms of Kronecker symbols, and one 
obtains for the local energy functional 



h{ai,<Tj;a) = -KS aiia . + L(S aua + 6 aj , a ) 

-MSa-aSv-a + E 2 . 



(12) 
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Collecting results the total hamiltonian reads 



the triangular lattice is given by |27ll3 



3 

The second line of the previous relation evaluates to a 
constant once the sum over all pairs and directions is 
performed, E tt E(ij)J £ (*<r ( 

) + E 2 } = N{E 2 + 

2E 3 ). The phase behavior is unaffected by this global 
shift of energy, thus, we discard the constant, and after 
rearranging terms, one arrives at the final expression for 
colloidal dimcr hamiltonian 



(ij) a=1 (ij)a 



(14) 



The excitation spectrum of orientational flips is governed 
solely by the two energy scales K, M to be compared to 
the thermal one ksT, and the phase diagram reduces to 
a two-dimensional plane. The K-term accounts for the 
energy gained if two neighboring dimers are in the same 
orientational state, see Fig.|5£b), £4 or E\. The M-term 
distinguishes the high symmetry state of alignment of the 
dimers with the lattice bond vector, Fig. |Hfb), E4. This 
coupling between spin configuration and orientation of 
the bond vector has a close analogue in the magnetic 
dipole-dipole interaction of localized spins on a lattice. 

The spin hamiltonian of Eq. Ijl4|l can be interpreted 
as a generalization of a g-state Potts model, the natu- 
ral extension of the Ising model. In these systems, each 
"spin" i can be in one of the o~i — 1, .., q equivalent states. 
The only energetic contribution occurs if two neighboring 
spins are in the same state, i.e., 



H 



(ij) 



(15) 



It is obvious that the colloidal dimer hamiltonian reduces 
to a three-states Potts model if M = 0. 

There have been intensive studies of the critical 
properties of ferromagnetic, J > 0, as well as anti- 
ferromagnetic, J < 0, Potts models for various lattice 
geometries in two and three dimensions. In particular, a 
rigorous solution of the ferromagnetic Potts model in two 
dimensions is available [29L |30( for a square, triangular, 
and honeycomb lattice for q — 2 (Ising) and q > 4. The 
melting of the ferromagnetic phase to the paramagnetic 
one occurs via a continuous phase transition for the case 
of q = 2,4 and via a discontinuous one for q > 4. For 
both cases the exact value of the critical temperature for 



exp^ = 2 cos (| cos- )) ; , = 2,4(16) 
eXP A^ = 2coshgcosh-(^)); 9 >4. 

Unfortunately, for the interesting case of q = 3 neither 
the method of transfer matrices (for the Ising model) nor 
the circle theorem |27j for the vertex model (for q > 4 
Potts models) is applicable. Nevertheless, since the crit- 
ical point in Eq. i fTfjj l agrees with the exact Ising result, 
it is expected that it holds also for q = 3. There is 
also strong numerical evidence by Monte Carlo simula- 
tions |30ll31| that the critical temperature is still given by 
Eq. (|16|) . Furthermore, Monte Carlo simulations as well 
as renormalization group analysis corroborate a continu- 
ous transition scenario [30( for q = 3. 

It is interesting to note that a standard mean-field 
analysis gives a qualitatively wrong result as far as the 
nature of the transition is concerned, i.e., a discontinuous 
scenario is predicted for q > 3. The value of the critical 
temperature is overestimated, 



3J 

2 In 2' 



(17) 



with respect to the tentatively exact value ksT c = 
J I ln[2 cos(7r/9)] by a factor of 1.3654 similar to the case 
of the 2D Ising model. Since we recover the three-state 
Potts model as a special case of the colloidal dimer hamil- 
tonian, Eq. HI 40 ! one expects that the mean-field analysis 
fails to predict the correct nature of the transitions also 
in the case of the colloidal dimer hamiltonian. 

For negative exchange coupling J with a low- 
temperature anti-ferromagnetic state there are no exact 
solutions known. For triangular lattices renormalization 
group (RG) studies [32, series analysis 0, and Monte 
Carlo simulations [34| give for the critical point a nu- 
merical estimate ksT c « J/ In (0.206). On the other 
hand, there is no general consensus on the nature of the 
transition. Extrapolated series expansions and Monte 
Carlo simulations find the transition to be discontinuous, 
whereas renormalization group analysis yields a continu- 
ous phase change |30l |. 



B. Mean- field description 

A rigorous solution of the colloidal dimer problem is 
certainly desirable, however, it would imply an analytic 
solution of the three-state Potts model as a special case, 
a problem that remains unsolved despite considerable ef- 
forts [271 |30j . Approximate methods are required to gain 
further insight in the the phase behavior. Mean-field the- 
ory has proven a powerful tool that allows to identify the 
global topology of the phase diagram, in particular, the 
coexistence of phases of different broken symmetries, as 
well as to provide approximate numbers for the location 
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of phase boundaries. The results of such a mean-field 
analysis also provide a useful reference for a numerical 
simulation by Monte Carlo methods which will be pre- 
sented in Sec. IIII CI 

We use the framework of variational mean-field theory 
and determine the optimal distribution of dimer states 
within a restricted class of density matrices Q . The nor- 
malized equilibrium density matrix p = Z^ 1 exp(— /3H), 
with inverse temperature f3 — l/k B T, is difficult to ob- 
tain since the partition sum Z = Tr [exp(— [3TL)] usually 
cannot be evaluated exactly. To overcome this problem it 
is favorable to approximate p by a product of single-site 
density matrices pi, 



p(Wi}) = 



(18) 



i.e., spins at different sites are taken as uncorrelated. 
The quantities Pi(<Ji) represent the probability to find 
the dimer at site i in state <ii. Normalization implies in 
particular Pii^i) = 1- The free energy correspond- 
ing to a system represented by a density matrix as in 
Eq. I|18|) results from the usual balance energy vs. en- 
tropy T p = Z p — TS P , where the mean energy £ p and 
the entropy S p are evaluated by the rules of statistical 
mechanics 



-fee (In p) p . 



(19) 



where (. . .) p denotes averaging with respect to the ap- 
proximate distribution of Eq. (|18|l . One can show in gen- 
eral that Bogolyubov's inequality Q sets a lower bound 
to F p , 



T p = Tr [pH] + k B TTr [plnp] > T . 



(20) 



Here T denotes the exact thermodynamic free energy. 
The idea of the variational mean-field method is to find 
the single-particle piipi) that yields the variational mean- 
field free energy T p closest to the thermodynamic one, 
viz. T . 

For the colloidal dimer hamiltonian the variational 
mean-field free energy reads 

T p = -k e E Pi(°)pj ( a ) - M E E p^ a )pi (°) 

(ij) <J a (ij) a 

+k B Tj2Y,p^ ln ^ a ^- ( 21 ) 



In the simplest case, all lattice sites are equivalent and 
the only phase transition allowed separates a paramag- 
netic phase from a ferromagnetic one. The suitable den- 
sity matrix is therefore sitc-indcpcndcnt, 



Pifa) = p{o-i) h-> Pa z 



(22) 



Upon substituting this ansatz in the mean-field free en- 
ergy, Eq. (|21[l. yields the free energy density, i.e., free 
energy per lattice site, / = J- p /N, 



f = -(3K + M) J2pI + k B T E P* ^ Pa- 



(23) 



Let us note that within the single-site mean-field descrip- 
tion, the parameters K and M appear only in the com- 
bination 3K + M. Applying the same procedure for the 
three-state Potts model gives an identical expression for 
the free energy density provided 3K + M is substituted 
by J. 

A ferromagnetic phase (FM) is characterized by a bro- 
ken symmetry, i.e., a particular dimer state is preferred 
with respect to the remaining. The single dimer exci- 
tations with respect to the ferromagnetic ground state 
are degenerate, i.e., the two possible minority orienta- 
tions are equally probable. By the symmetry of the FM, 
this property is preserved in the whole phase region. Se- 
lecting the majority component as the a = 1 state, we 
parametrize p\ = (l + 2m)/3 by a generalized magnetiza- 
tion m. The minority probabilities are determined by the 
normalization J2 a P<? ~ t° P2 = P3 = (1 — ro)/3. The 
magnetization is restricted to —1/2 < m < 1, and the 
fully aligned ferromagnetic state corresponds to m = 1, 
whereas the paramagnetic one is given by m = 0. A 
negative magnetization would indicate that the distin- 
guished spin state is less probable than the remaining 
two. Hence there is an inherent asymmetry m i— *■ — m in 
the mean- field description, and as we shall see later, all 
transitions within this mean-field approach are (falsely) 
predicted to be discontinuous, similar to the case of a 
pure Potts model. 

Substitution into Eq. (|23|l yields the excess free energy 
of the ferromagnetic phase with respect to paramagnetic 
one, 



/fm - fp 
k B T 



2c 



1 + 2m 



ln(l + 2m) 



„ 1 — m , , , 
+2— — ln(l-m) ; 



(24) 



Note that in the M-term the spin orientation and the lat- 
tice bond are aligned, a = a. This expression for the free 
energy is still subject to the minimization procedure with 
respect to appropriate sets of single-site density matrices. 

The symmetries of the dimer states in the broken sym- 
metry phase have to be reflected in the ansatz for the 
Pi{o~i). In general, one introduces a number of sublat- 
tices each comprising a set of equivalent lattice sites with 
respect to orientational order. 



which depends on the effective parameter e = f3(3K+M). 
The corresponding graph is shown in Fig. |}5j(a). One in- 
fers the typical characteristics of a first order phase tran- 
sition, e.g., a discontinuous change of the order parameter 
and metastable states. 

The equation of state for the magnetization is obtained 
by minimizing /fm with respect to m, 



2em = In 



1 + 2m 
1 — m 



(25) 
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T<T C 



T< T c 




m 

FIG. 9: a) Free energy as a function of the magnetization for 
various temperatures above, below, and at the critical point. 
Local minima correspond to metastable phases, b) Mean-field 
temperature dependence of the magnetization for the ferro- 
magnetic phase of the colloidal dimer hamiltonian. At the 
critical temperature, T/T c = 1, the order parameter exhibits 
a jump. 



One infers that m = 0, corresponding to the paramag- 
netic phase, is always a solution to this equation. De- 
pending on the value of the effective parameter e non- 
trivial solutions may occur that exhibit lower mean- 
field free energies than the paramagnetic one, /p = 
—ksT(e/3 — In 3). For infinite temperatures (e = 0) 
the solution to this equation reduces to the paramagnet 
(to = 0), whereas at zero temperature (e — > oo) full align- 
ment of the ferromagnet (to — 1) occurs. The complete 
temperature dependence of the magnetization is depicted 
in Fig. El (b). 

The thermodynamic transition from a paramagnetic 
phase occurs at the temperature T c when the non-trivial 
minimum of the free energy /fm equals the one of the 
paramagnetic state fp. With the help of the equation of 
state, Eq. (|25Jl , one finds for the critical effective parame- 
ter e c = 2 In 2, which implies for the critical temperature 



FM - P 



3K + M 
21n2 



(26) 



Representing the phase diagram in terms of the dimen- 
sionless parameters k B T/\K\ and M/\K\, the FM-P 
phase boundary yields a straight line, see Fig. ^3 

The magnetization at the critical point is finite and 
attains the value m c = 1/2. The jump in the derivative 
of free energy density with respect to temperature at the 
phase transition implies a latent heat per dimer which 
evaluates to qi/ksT c = ln(2)/3. The heat capacity jumps 
at T c by Ac ^ = (8fc B /3)[ln^ (2)/(3 -41n(2))]. Note that 
these quantities are universal in the sense that they do 
not depend on the parameters K and M, i.e, they are 
constant along the whole phase boundary. The mean- 
field analysis appears to be qualitatively correct for K < 
0, whereas for K > the simulation results suggest that 
fluctuations render the phase transition to a continuous 



one. as will be discussed in detail in Sec. IIII CI 
K>0 K<0 




k B T/K 



FIG. 10: Phase diagram for colloidal dimers. Solid lines rep- 
resent first order transitions derived from the MF description 
and the dashed line after the tricritical point represents a line 
of second order transitions. Dots and triangles denote the 
respective continuous and discontinuous transition points as 
obtained by MC simulations; the two encircled symbols in- 
dicate the pure Potts transitions. Dotted lines represent the 
extrapolation of the critical point of the Potts model, see text. 
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FIG. 11: Ground state ordered structures of 2D colloidal 
dimers: (a) ferromagnetic (FM), (b) herring bone (HB), 
(c) Potts anti- ferromagnetic (AM), and (d) Japanese 6 in 1 
(J6/1), structure. In the latter, the energy of the system does 
not depend on the orientation of dimers on sites denoted by 
circles. 

The ferromagnetic phase represents only the simplest 
broken symmetry. Spatially varying order parameters 
give rise to complex phase behavior, in particular, the 
sites fall into different classes of sublattices. Compar- 
ing the exact ground state energies of different ordered 
structures, one can infer the possible broken symmetry 
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phase at low temperatures. In total, we have found four 
ordered phases of different symmetry, see Fig. ^2 that 
are thermodynamically stable in different regions of the 
phase diagram, depicted in Fig. In addition to the 
ferromagnetic one, there appears an anti-ferromagnetic 
phase (AM) stable for K < and not too large M. Here, 
in the fully ordered state, neighboring dimers are in dif- 
ferent orientational states. On the triangular lattice this 
implies that each dimer is surrounded by three neighbors 
of a configuration rotated by 60° and by three neighbors 
of a configuration rotated by —60° with respect to the 
configuration of the central dimer. There are three tri- 
angular sublattices with a lattice constant of aV3 and 
the sublattices are equivalent in the sense that each of 
them can be obtained by a spatial translation with a si- 
multaneous rotation of the dimers of any other one. For 
positive K and sufficiently negative M a herring bone 
(HB) structure represents the state of lowest free energy, 
see Fig. ^2 Statistically equivalent dimers are arranged 
in two sets of rows, the majority orientational states cor- 
responding to the ones not parallel to the orientation of 
the rows. Finally, we have found that intervening in be- 
tween the FM and AM for K < there appears a phase of 
four triangular sublattices. Three of them are equivalent 
in the sense mentioned above, the fourth one exhibiting 
no preferred orientation. We refer to it as Japanese 6 
in 1 phase (J6/1), due to its resemblance to a weaving 
pattern for chain mailles worn by samurais in the 14 th 
century p35| . 

To describe these complex phases generalized ansatzes 
have to be used for the density matrix. Then the site 
dependence of the components of pi(oi) is encoded in 
the sublattice T the site i belongs to, i.e., 



Pi{cri) 



i € sublattice T . 



(27) 



Some of the density matrices p\ are related by permu- 
tations of dimer states. Explicitly, for the AM phase as 
well as for the HB the sublattices are equivalent, so that 
only one density matrix characterizes the order of the sys- 
tem. For the J6/1 phase the order in the three equivalent 
sublattices is determined by one non-trivial density ma- 
trix, whereas the fourth one attributes equal probabilities 
to all dimer orientations. Let us mention an additional 
symmetry that occurs for the AM and J6/1 phase. Here, 
the two minority orientations are equivalent as was the 
case in the FM phase. The corresponding density ma- 
trix is parameterized by a single number, a generalized 
magnetization m defined on a sublattice. The situation 
is different for the HB phase. An excitation of a dimer 
to the majority orientation of the neighboring row is in 
general different from aligning the dimer with the orien- 
tation of the row. Correspondingly, the density matrix 
in the HB phase requires two parameters. 

The three sublattices for the anti-ferromagnetic phase 
are inferred from the ground state represented in Fig. 1111 
The single-site density matrices corresponding to the sub- 



lattices A, B,C, see Eq. (|27() . are parameterized by 
P 2 B =P 3 C -(l + 2m)/3, 



pi 

P2,3 



(28) 



The excess free energy density of the anti-ferromagnet is 
then obtained by substituting this ansatz into Eq. (|21|) . 



/am — / ] 



p 2u 9 1 + 2to , . 

= — h-m 2 + —, ln(l + 2m) 



+2- 



3 3 
1 — m 



In (1 — m) 



(29) 



where the sole parameter entering this expression is given 
by (i = —0(3K + M)/2. The shape of the free energy 
density /am as a function of temperature is identical 
to the one of the ferromagnetic phase /fm, see Fig. El 
The equation of state for the sublattice magnetization 
to is thus obtained by replacing e by p in Eq. I|25|). i.e., 
2pm = ln[(l + 2m)/(l — to)]. The same substitution for 
the critical temperature, Eq. i|26|) . yields 



AM-P 



k B T r 



3K + M 
41n2 



(30) 



The AM-P phase boundary is represented by a straight 
line in the T/\K\—M/\K\ phase diagram, Fig EH The 
latent heat and the excess of the heat capacity are again 
the same as for the FM phase provided one substitutes e 
with p. 

For the Japanese 6 in 1 phase the appropriate ansatz 
for the three equivalent sublattices (A,B,C) is identi- 
cal to the anti- ferromagnetic case, Eq. (|28() . whereas the 
fourth (O) is degenerate 



P°2,3 = l/3- 



(31) 



The excess free energy density for the J6/1 symmetry is 
readily obtained 



/. 



J6/1 



fp 



1 + 2to 



-TO 



3 

1 — TO 



In (1 + 2m) 



In (1 — to) 



(32) 



Here v = —f3(3K — M)/3 denotes the sole effective pa- 
rameter. Up to the multiplicative constant 3/4 the ex- 
cess free energy density of the J6/1 phase is identical 
to the excess free energy densities of the FM and AM 
phases discussed above, provided v is substituted by e or 
p, respectively. The prefactor 3/4 arises since only three 
of the four sublattices exhibit a broken symmetry. The 
equation of state obtained as the derivative of the free 
energy density with respect to the sublattice magnetiza- 
tion to is again formally identical to the cases discussed 
above, 2^to = ln[(l + 2m)/(l — to)]. The critical temper- 
ature for the J6/1-P phase transition is readily obtained 
via the substitution rule, yielding v c = 2 In 2 and 



J6/1 



knT c 



3K-M 
61n2 



(33) 
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The J6/1-P phase boundary is indicated in Fig.^H Since 
only three of the four sublattices exhibit a broken sym- 
metry, the latent heat as well as the jump of the heat 
capacity at the phase transition are 3/4 of the universal 
values of the FM-P or AM-P transitions. 

The herring bone structure (HB) requires two parame- 
ters, one for the majority magnetization to and a second 
'biaxiality' parameter, n, for the asymmetry of the minor- 
ity orientations. The need for the biaxiality parameter 
becomes clear if one considers excitations of HB ground 
state. A single dimer flip to the orientation of the neigh- 
boring row costs an energy of AE ~ —2K — 2M, whereas 
aligning it with the row induces an energy change of only 
AE = 2K. (Note that the single dimer excitation ener- 
gies become degenerate for M = — 2K.) The low temper- 
ature excitations of the HB structure are represented by 
a dilute gas of defects of the latter kind. Considering the 
sign of the single-dimer excitation energies an estimate 
for the stability of the HB structure at low temperatures 
is possible, i.e., for positive K and sufficiently large neg- 
ative M. 

For the parameterization of the sublattice single-site 
density matrices we choose 



92 = Pi = (1 + 2m)/3 , 
Pi,3 = Pi,2 = (1 - m ± n)/3 . 



(34) 



Then the free energy density corresponding to the HB 
symmetry reads 



/hb — ,/p 



k B T 3k B T 



M 2 «|I , 
-to n 



9k B T 



2(2K + M) 



l + 2m 





3k E 


T 


1 


— to 


+ B 1 




3 




1 


— to 




3 



In (1 — m + n) 
In (1 — m — n) . 



In (1 + 2m) 



(35) 



The excess free energy is subject to a simultaneous min- 
imization with respect to both order parameters to, n. 
In general, one cannot find an analytic solution of the 
equations of state 



1 — if) — I— n 2 

k B Tbx — = -2C2K + M)m+ -{AK + M)n, 

1 — to — n 3 



k B T\n 



(1 + 2to) 2 
(1 — to) 2 — n' 2 



= -2Mm - 2{2K + M)n . (36) 



The generic behavior close to the HB-P phase transition 
is illustrated in Fig. El f° r a path of constant M/\K\ 
close to the HB-FM line and one starting deep in the 
HB phase. 

For the case of a continuous transition the transition 
temperature can be obtained analytically by expanding 
the excess free energy density in powers of the small quan- 




FIG. 12: Order parameters m and n as a function of tem- 
perature for a path of constant M/K. (a) A discontinu- 
ous phase transition characterized by a jump of m and n 
for M/K = —2.1 and (b) a continuous phase transition for 
M/K = -5. 

tities to and n, which up to quadratic order reads 
/hb-/p f M An 4K + 2M 



k B T 



+ 1 TO + 

3k B T J 3k B T 
4JC + M\ n 2 
3k B T ) Y ' 



1 



(37) 



The quadratic form ceases to be positive definite and the 
paramagnetic phase becomes unstable via a continuous 
transition at the critical temperature 



HB - P, M < 



IK 



kuT r 



2{K + M) 



(38) 



and the critical direction is along the diagonal (to, n) ~ 
(1,1), i.e the order parameters become equal close to 
the transition, see Fig. In the regime —2K < 

M < —7K/3 the melting of the HB structure towards 
the paramagnetic phase is discontinuous. The tricritical 
point, where the order of the phase transition changes 
can by determined again by local considerations. A 
suitable non-linear transformation of variables (to, n), 
shows that the origin becomes unstable with respect to 
higher order terms along the critical direction at M = 
-7K/3, k B T 3c = 8K/9. 

For M = —2K the excitation energies of the isolated 
single-dimer flips in the HB structure become degener- 
ate. This is reflected by the fact that the excess free en- 
ergy, Eq. (|35|l . exhibits the additional symmetry n — > — n 
along this line. At zero temperature the HB phase is 
perfectly ordered (to = \,n — 0), i.e in a n-symmetric 
phase. Increasing the temperature does not give rise to 
a continuous or discontinuous n-symmetry breaking, so 
n = holds along the M = —2K line for all tempera- 
tures. With this observation in mind, one checks that 
the mean-field free energy density of the HB phase for 
M = —2K equals the one of the FM phase for all temper- 
atures. Thus, the line M = —2K is a coexistence line of 
the two neighboring phases. The special point for the dis- 
continuous transition to the P phase, k B T c = K/(2\n2), 
constitutes a triple point where FM, HB, and P phases 
coexist. 

The remaining coexistence lines between different or- 
dered phases, i.e., for K < between FM and J6/1 and 
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TABLE I: Phase boundaries obtained in mean-field theory. 
The approximate results for the J6/1-FM and J6/1-AM tran- 
sition correspond to the asymptotic low-temperature (high- 
temperature) solution. 



TABLE II: Mean-field special points. 



transition 


critical line 


K > 




FM-P 
HB-P 
r M— rirJ 


k B T c = (3Jsr + M)/(21n2) 

k B T c = -(2/3) + M), M < -7K/3 

M — —IK 


K < 




FM-P 
J6/1-P 
AM-P 
J6/1-FM 

J6/1-AM 


k B T c = (3K + M)/(21n2) 

k B T c = (-3K + M)/(61n2) 

k B T c = -(3K + M)/(41n2) 

k B T c » 3(5/sT + M)/(21n2), M < 6\K\ 

k B T c » 2(5if + M)/ln3, M > 5|K| 

fc B T c » -3(tf + M)/(21n2), M > 3\K\/5 

k B T c w -2(if + M )/ In 3, M < | Jf| 



between AM and J6/1, can only be determined numer- 
ically and are included in the phase diagram. Here we 
present in short approximate analytical solutions in terms 
of asymptotic low-temperature expansions. 

At the coexistence line the free energy densities of the 
two coexisting phases are equal. In particular, for the 
FM-J6/1 line /fm ( m FM ? T) = /j6/i(m J6 / 1} T), where 
/fm and /j6/i are defined in Eqs. ((24) and 1 (52) 1 and the 
m's are solutions of the corresponding equations of state. 
For T -> (M/|if| > 5), m FM ,mj 6 /i « 1 which yields 
the coexistence line 



FM - J6/1, M/\K\ > 5 : k B T c 



2(5K + M) 
in3 



(39) 



Close to the FM-J6/1-P triple point, M = 6\K\, 
k B T = 3|if |/(21n2), the magnetizations approximately 
take their critical values, n r iFM,n r ij6/i ~ "ic = 1/2, and 
give rise to 

FM-J6/1, M/\K\<6: k B T c = ■ (40) 



For the AM-J6/1 line an equivalent procedure with the 
same approximate solutions for the magnetizations close 
to the ground state or triple point results in 

AM - J6/1, M/\K\ < 1 : k B T c a ; (41) 



Analytical expressions for exact or approximate critical 
lines are gathered in Tab. [fl and the mean-field special 
lines, i.e., tricritical and triple points, in Tab. ITTl 



coexisting phases 


[k B T/\K\,M/\K\] 


K > 




FM-HB-P 

HB-P tricritical point 


[1/(2 In 2) -2] 
[8/9,-7/3] 


K < 




FM-J6/1-P 
AM-J6/1-P 


[3/(2 In 2), 6] 
[3/(5 In 2), 3/5] 



C. Results of the Monte Carlo simulations 

In this subsection we exemplify the phase behavior of 
the colloidal dimer system by numerical exact means, i.e., 
Monte Carlo methods. In particular, we corroborate the 
existence of all the phases discussed within the mean-field 
analysis and address the nature of the phase boundaries 
separating the corresponding phases of broken symmetry. 
The results obtained within mean-field thereby play a 
crucial role to identify the relevant order parameters to 
be monitored. 

The standard Metropolis algorithm [3j| has been em- 
ployed to generate a sequence of configurations repre- 
sentative for the canonical ensemble to some specified 
temperature. The starting configuration for the lowest 
temperature has been chosen as the perfectly ordered 
(exact) ground state corresponding to the pair of en- 
ergy scales (K,M). The temperature has been increased 
gradually starting from the representative configuration 
of the previous parameter set. We have considered sys- 
tems of different sizes subject to periodic boundary con- 
ditions. Explicitly, the number of rows L, each consisting 
of L dimcrs, has been varied, L = 12, 32, 52, 102, to en- 
sure that the results have approached the infinite system 
size limit. Before actual data were collected the system 
has been carefully equilibrated by performing up to sev- 
eral thousand Monte Carlo cycles. Each Monte Carlo cy- 
cle consists of NlnN attempted dimer flips, where 
N = L 2 is the total number of dimers in the system. To 
ensure that the configurations entering the data analysis 
are statistically independent, we have considered config- 
urations separated by up to 1000 Monte Carlo cycles in 
the vicinity of a phase transition. 

To characterize the broken symmetry phases we have 
considered the respective order parameters. For example, 
in the ferromagnetic case the fluctuating magnetization 
reads 



M = — - 1 ^ a ^ s 
N 2 N 



(43) 



Here, M a is the fluctuating number of dimers in the state 
(7. In particular, A/e corresponds to the majority orien- 
tation in the given configuration, and N — Yl^N'g is the 
total number of the dimers in the system. The macro- 
scopic magnetization m = (A4) is evaluated by averaging 
over the statistically independent configurations obtained 
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as discussed above. We have also measured the variance, 
((5M) 2 ), 5M = M. — m, which is related to a generalized 
susceptibility 

X=T^((M-mf), (44) 
k B T 

by a fluctuation-dissipation theorem. Furthermore, we 
have computed the average energy U = (H) (with u — 
U/N we will denote the average energy per dimer) as well 
as the corresponding variance related to the specific heat 
per dimer 

For the other broken symmetries the fluctuating sublat- 
tice magnetizations are defined to reproduce the defini- 
tions given in the previous Section. 

In the following, the phase behavior as well as the 
phase transitions for the different cases will be illustrated. 
Since a mean-field result does not always predict the cor- 
rect qualitative behavior, in particular the order of the 
phase transition, see Section [ill Bl we have put consider- 
able effort to investigate the nature of the phase bound- 
aries. 

1, Ferromagnetic system 

By comparison of the free energies of different ordered 
configurations we have found that the ground state with 
the ferromagnetic broken symmetry is realized for K > 
with M/K > -2 and for K < with M/\K\ > 5. Due 
to the different mechanism of the ferromagnetic ordering 
in the two cases, the corresponding phase transitions to 
the paramagnetic phase differ in their nature and charac- 
teristics. In the following we will discuss the two regions 
separately. 



(a) (b) 
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FIG. 13: The ferromagnetic system (K > 0, M/K = 1). Tem- 
perature dependence (a) of the order parameter and its vari- 
ance and (b) of the energy per dimer and its variance. Both 
indicate a continuous nature of the phase transition. 

(i) For K > the parallel ordering of dimers is fa- 
vored by the K-tetm of the colloidal dimer hamiltonian, 



Eq. i|14|) . i.e., by its Potts part. In addition, the remain- 
ing M-term also favors this preference and effectively 
rescales the energy scale K in the same way as it has 
resulted from the mean-field description. In essence, in 
this regime the colloidal dimer system can be described 
by an effective hamiltonian, 7i c s = — Kgg <Wo-i) 
where K e s = K + M/3. Our Monte Carlo simulations 
are in agreement with the known results for such a hamil- 
tonian. The FM-P phase transition is clearly continuous 
as can be inferred from the pronounced fluctuations of 
the magnetization and energy in the vicinity of the tran- 
sition, see Fig. El Furthermore, the distribution func- 
tions of the magnetization and energy per dimer exhibit 
only a single peak in the whole temperature range under 
consideration (not shown). 

The critical temperature of the FM-P transition 
closely follows the "exact" value of the Potts hamiltonian, 
Eq. (|16f> , provided the exchange coupling J is replaced by 
the above mentioned effective energy scale, K e g. Only 
close to the coexistence region with the HB phase, for 
— 1 > M/\K\ > —2, the critical temperature deviates 
from the scaling law in K c g and the ferromagnetic order 
melts before the Potts critical temperature is reached, 
see Fig.llUI Nevertheless, within the numerical precision 
the nature of the transition remains unaffected. 

(ii) For K < the Potts term of the colloidal dimer 
hamiltonian, Eq. Q14[). disfavors ferromagnetic order and 
anti-ferromagnetic order would be expected. However, 
for high values of the energy scale M, explicitly for 
M/\K\ > 5, the term favoring the alignment of neigh- 
boring dimers along the line connecting them dominates, 
and the system exhibits a FM symmetry at low temper- 
atures. 



(a) (b) 
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FIG. 14: The ferromagnetic system (K < 0,M/\K\ = 5.5). 
Temperature dependence (a) of the order parameter and its 
variance and (b) of the energy per dimer and its variance 
(right). Both indicate a discontinuous nature of the phase 
transition. 

The critical properties as well as the mechanism of the 
transition to the P phase are expected to differ from those 
of the Potts hamiltonian. Indeed, the MC simulations of 
the full hamiltonian reveal that the phase change is dis- 
continuous. At the critical point the systems "jumps" be- 
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tween the two degenerate, FM and P, phases, and typical 
snapshots reveal phase coexistence regions at the critical 
point (not shown) . The discontinuous nature of the tran- 
sition is also evidenced by the fluctuations of the order 
parameter and energy, see Fig. 1141 

The occurrence of metastable phases can be quantified 
by determining the distribution functions for the fluc- 
tuating values of the order parameter and energy per 
dimer [37|. We have carefully collected their distribu- 
tions in the vicinity of the transition and confirm a dis- 
continuous transition. 



2. Herring bone phase 

The broken symmetry state with a herring bone pat- 
tern is realized for K > and M/K < -2. Within 
mean-field theory we have shown that its ordering is de- 
scribed by two order parameters, i.e., the magnetization 
to and the biaxiality parameter n. In MC simulations we 
have monitored them both. 
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FIG. 15: The herring bone system (K > 0,M/K = -2.1). 
Temperature dependence (a) of the order parameter and its 
variance and (b) of the energy per dimer and its variance. 
The distribution functions (c) of the order parameter and (d) 
energy per dimer. At the phase transition, note the bimodal 
structure characteristic for a discontinuous phase transition. 

Qualitatively, the MC temperature dependences of m 
and n are the same as the ones determined within the 
mean-field description. We find a discontinuous phase 
change for M/K < —2, identified by jumps of the mag- 
netization and energy per dimer at the critical temper- 
ature. The jump of the biaxiality parameter is masked 
by rounding due to the finite size of the system. How- 



0.1 - 



03 
Q. 



k B T/K= 1.00 
]T/K = 0.58 
k B T/K = 0.34 




-10 



-6 
M/K 



FIG. 16: The probability of finding an M-pair [p(Af-pairs) = 
Nm -pairs /3N] in a herring bone system in the HB phase 
(ksT/K = 0.34), close to the saturated value of the critical 
temperature (kBT/K = 0.58 and ksT c /K ~ 0.6), and well in 
the P phase (k B T/K = 1) as a function of M/K. Note that 
the equilibrium probability for M-pairs in the normal P phase 
(where the probabilities of all dimer orientations are equal, as 
well as the probabilities of all configurations of dimer pairs) 
is 1/9. 



ever, we find a sharp increase upon approaching the crit- 
ical point from below, similar to the one found within 
the mean-field description, see Fig. 1151 For increasing 
negative values of M/K the latent heat decreases, so 
that below M/K = —4 we cannot distinguish between a 
weakly discontinuous and a continuous transitions. Cor- 
respondingly, it is difficult to clearly corroborate the ex- 
istence of the tricritical point predicted by mean-field at 
M/K = —7/3, let alone to determine its position. Mon- 
itoring the probability distribution of the order parame- 
ters suggests that indeed there is a change of the nature 
of the phase transition in the vicinity of the mean-field 
prediction. 

It is known that quantitatively, the effect of fluctua- 
tions decreases the critical temperature obtained within 
mean-field. As it can be clearly seen in Fig. ^3 for the 
HB-P transition this appears to be not the only effect. 
If in the interval — 2 > M/K > — 3 there is still some ev- 
idence of the scaling law similar to the one found within 
the mean-field description, this practically vanishes for 
M/K < -3. Finally, for M/K < -4 the critical temper- 
ature "saturates" at a value ksTc/K « 0.6, i.e., it does 
not depend any more on the energy ratio M/K. 

To rationalize this finding, let us recall the expression 
for the colloidal dimer hamiltonian l|14|) . The first, K- 
tcrm, favors for K > parallel alignment of dimers. On 
the other hand, the M-term with a large negative value 
of M strongly suppresses parallel alignment in which the 
two dimers face each other - to be referred to as an M- 
pair, Fig. 02b). As a consequence, the probability of 
finding M -pairs in either of the two phases, HB or P, 
is rapidly decreasing with increasing negative values of 
M/K. The number of M-pairs as a function of M/K 
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is shown in Fig. ^] for three paths of a constant tem- 
perature, viz. in the HB phase, in the P phase, and at 
the intermediate temperature ksT/K = 0.58 for which 
the system traverses the phase boundary from HB to P 
at M/K » —3.7. For large enough negative M the role 
of the M-term is only to reduce the number of degrees 
of freedom in the system. Then the corresponding criti- 
cal behavior depends solely on the relative magnitude of 
the exchange energy K with respect to the thermal en- 
ergy ksT. In the {UbT/\K \,M/\K\) phase diagram such 
a phase boundary is represented by a straight vertical 
line. Within mean-field the reduction of the number of 
degrees of freedom for the configurations of dimer pairs 
is not properly captured. Thus the role of the M-term is 
as in all the other cases reduced to effectively rescale the 
exchange energy K . 



3. Anti-ferromagnetic phase 

The ground state of the anti-ferromagnetically ordered 
phase occurs in colloidal dimer systems with K < 0, 
M/\K\ < 1. As discussed in Sec. IIII Bl its order can 
be described by a single order parameter - a sublattice 
magnetization m, where the system is divided into three 
triangular sublattices. 
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FIG. 17: The anti-ferromagnetic system (K < 0,M/\K\ = 
0.4). Temperature dependence (a) of the order parameter and 
its variance and (b) of the energy per dimer and its variance. 

The AM-P phase transition is discontinuous in the 
whole range of the stable AM phase. In contradiction to 
the universal jump of the latent heat and heat capacity 
determined by the mean-field approach, the MC results 
show slightly decreasing strength of the discontinuity of 
the transition when M/\K\ is increased towards 1. At the 
special point, M = 0, where the colloidal dimer hamilto- 
nian reduces to the Potts model, our results match the 
ones known for the Potts anti-ferromagnet [30l l34j . Ap- 
plying again the empirical shift of the mean-field analysis 
to match the anti- ferromagnetic Potts point, gives a satis- 
factory description of the phase boundary for M/\K\ > 1, 
see Fig. 1101 For highly negative M/\K\ the critical tem- 
perature saturates as is the case of the HB-P transition. 



The saturation again originates in the suppression of M- 
pairs as already discussed above. 



4- Japanese 6 in 1 system 

For K < 0, squeezed between the FM and AM phases 
one observes the new J6/1 phase. Its highly degenerate 
ground state is stable for 1 < M/\K\ < 5. The degener- 
acy of the J6/1 phase originates in the disordered nature 
of one of its sublattices. A typical snapshot of a ground 
state in the J6/1 phase is presented in Fig. In addi- 
tion, a schematic representation and a picture of a chain 
maille using the Japanese 6 in 1 weaving pattern |35| 
from which we have borrowed the name are plotted. 
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FIG. 18: The Japanese 6 in 1 system. Left: A snapshot of one 
of its ground states. Middle: Schematic representation of the 
J6/1 configuration. Here the circles indicate the degeneracy 
of the dimer orientation on the disordered sublattice. Right: 
A picture of chain maille using the Japanese 6 in 1 weaving 
pattern [BE!. 
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FIG. 19: The Japanese 6 in 1 system (K < 0,M/\K\ = 3.5). 
Temperature dependence (a) of the (sublattice) order param- 
eter and its variance and (b) of the energy per dimer and its 
variance (right). Both indicate a continuous nature of the 
phase transition. 

The MC simulations corroborate the division of the 
J6/1 system into four triangular sublattices. Three of 
them are equivalent and ordered with mutually different 
preferred dimer orientations for temperatures below the 
critical one. The remaining sublattice exhibits no pre- 
ferred orientation for all temperatures. The magnetiza- 
tion parameter is defined as the average of the three non- 
zero sublattice magnetizations, M. — ^ r A4r/3. In the 
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FIG. 20: The system (K < 0,M/\K\ = 5.2) exhibits a se- 
quence of phase transitions FM-J6/1-P upon increasing the 
temperature. The FM-J6/1 structural transition takes place 
at T c .i and is discontinuous whereas there is no jump in the 
order parameter and energy per dimer at the J6/1-P transi- 
tion, which takes place at T c ,2- 



simulation we have monitored both, the magnetization of 
the ordered as well disordered sublattice. Upon increas- 
ing the temperature from the ground state, the order in 
the ordered sublattices of the J6/1 phase melts and the 
system exhibits a continuous phase transition to the P 
phase. The corresponding phase boundary is represented 
by a curved, non-monotonic line, see Fig. EH Close to the 
critical point the fluctuations of the order parameter are 
strongly enhanced and the system can "switch" between 
the four classes of configurations, i.e., characterized by a 
translation of the disordered sublattice. Graphs of the 
typical temperature dependences of the sublattice mag- 
netization and energy per dimer are depicted in Fig. I19I 
The high residual entropy stabilizes the J6/1 with re- 
spect to the neighboring FM or AM phase upon increas- 
ing temperature. Correspondingly, the parameter region 
where the J6/1 phase is stable expands upon heating. 
For parameters 5 < M/\K\ < 5.3 or 0.8 < M/\K\ < 1 
the melting of the FM or AM proceeds in two steps via a 
discontinuous structural transition to the J6/1 phase fol- 
lowed by another transition to the paramagnetic state. 
A typical temperature variation of the energy per dimer 
and the corresponding heat capacity for the case of a FM- 
J6/1-P sequence is presented in Fig. EH The AM-J6/1- 
P phase sequence is characterized by similar graphs. 



IV. CONCLUSIONS 

We have shown that soft matter colloidal systems ex- 
hibit peculiar phase behavior once exposed to external 
laser fields. The laser-colloid interaction favors configu- 
rations where the particles are attracted to the spots of 
high light intensities. By a subtle balance of the com- 
pression of groups of two or three particles at such spots 
and their mutual repulsion composite objects are formed, 
such as dimers and trimers. The relevant low-energy de- 



grees of freedom are then the discrete orientational con- 
figurations with an excitation spectrum at the thermal 
scale. 

The introduction of an effective hamiltonian allows to 
map the phase behavior to spin-like systems. In partic- 
ular, we have rationalized the case of trimers on a tri- 
angular lattice [2^[. which has motivated this work, in 
terms of a two-dimensional Ising model. A comparison 
with experimental data corroborates the approach of us- 
ing composite objects and validates our theory. 

We have suggested a new experimental setup where 
dimers play the role of the composite objects. It is 
shown that a rich phase behavior emerges; in terms of a 
"spin language" a ferromagnetic, an anti-ferromagnetic, 
a herring-bone, and an exotic Japanese 6 in 1 phase have 
been determined. We have mapped the phase diagram 
using a variational mean-field theory and substantiated 
our findings by Monte Carlo simulations. 

It is interesting to note that these colloidal molecular 
crystals have a counterpart on the atomic scale, viz. ad- 
sorption of molecules on atomic surfaces. For instance, 
nitrogen on graphite |38t l39| at suitable densities orders 
in a herring-bone structure similar to the one of colloidal 
dimers, which lead Sluckin (Ijj to construct a general- 
ized Potts model closely along the line of reasoning pre- 
sented here. Allen and Armitstead [ilj presented Monte 
Carlo simulations close to the transition from ferromag- 
netically ordered nitrogen adsorbats to the experimen- 
tally observed herring-bone structure. Using solely sym- 
metry arguments, Vollmayr |42j discussed the possible 
ground states if the interaction between two molecules 
depends on the bond orientation and explored regions of 
suppressed relaxations in an a priori unfrustrated sys- 
tem. 

It appears promising to extend the concept of com- 
posite objects to different symmetries of the underlying 
lattice, e.g. square or rectangular unit cells, and to non- 
integer filling factors. The corresponding phase diagrams 
are expected to include a variety of phases and phase 
transition. 
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APPENDIX A: RELATION BETWEEN ENERGY 
SCALES AND EXPERIMENTAL CONTROL 
PARAMETERS 

In this Appendix we show how the parameters of the 
trimer hamiltonian can be expressed in terms of the laser 
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intensity, the screening length of the colloid-colloid inter- 
action, and the effective charge of the colloids. The case 
of dimers is then a straightforward modification and we 
merely state the result at the end of this Appendix. 

The formation of composite objects, i.e., trimers, as 
well as their mutual interaction is due to a subtle in- 
terplay of the external laser potential and the screened 
Coulomb colloid-colloid repulsion. The symmetry of 
the external potential favors equilateral triangles aligned 
along one of the two distinguished lattice directions. The 
potential squeezes the constituents of the trimer closer 
to each other. The forces of the external potential on a 
trimer are balanced by the internal constituent Coulomb 
repulsion. To lowest order the effect of the neighboring 
trimers can be ignored since in the regime of interest the 
trimers are small compared to the lattice constant. Fur- 
thermore the Coulomb interaction is exponentially sup- 
pressed for the case of screening lengths smaller than the 
lattice constants. The thermal pressure is negligible com- 
pared to the binding energy. 

To this order of approximation the ground state of the 
system has a two-fold degeneracy per lattice site. Intro- 
ducing the residual interaction of the trimers, the energy 
levels are split on a scale comparable with the thermal 
one. 

Let us discuss the interactions of the trimer con- 
stituents in detail. First the dielectric contrast at optical 
frequencies of the colloidal particles with respect to the 
solvent gives rise to a coupling to the laser intensity. This 
interaction can be expressed in terms of an effective po- 
tential. For the experimental setup under consideration 
the total electric field consists of a coherent superposition 
of three laser beams 

E(r, t) = Eoe - *"* [e lkl ' r + e jk2 ' r + e lk3 ' r ] + c.c, (Al) 

where the respective in-plane projections of wave vectors 
k constitute an equilateral triangle. The total intensity 
consists of a rapidly varying component in time, which 
averages to zero. The remaining time-independent com- 
ponent gives rise to a constant, which we also discard, 
and the interesting part, characterized by a spatial mod- 
ulation, 

7(r) = 4|£ | 2 [cos(Gi ■ r) + cos(G 2 ■ r) + cos(G 3 • r)JA2) 

where Gi = k 2 — k 3 , G 2 = k 3 — ki, and G3 = ki — k 2 
again form an equilateral triangle. 

The periodic potential exhibits the symmetry of a tri- 
angular lattice with G»'s being the wave vectors of the 
reciprocal lattice. In the following we choose the coordi- 
nate system in which 

G 1 = G(0,-1), G 2 = G(^,i), G 3 = G(-^,i). 

" . (A3) 

It is also convenient to introduce the spacial lattice vec- 
tors 

ai=o(l,0), a 2 = a(i,^), as = a{-\, % .(A4) 



where a = 4n/(Gy/3) is the lattice constant of the exter- 
nal potential. In this parameterization and with respect 
to the lattice sites the relative positions of the colloidal 
particles within a trimer with orientation S G {+1,-1} 
are parameterized by the single dimensionless distance 
parameter p: 

Ri = Spai, R 2 = S>a 2 , R 3 = S pa 3 . (A5) 

The effective interaction of the colloidal particles with 
the external electric potential is proportional to the laser 
intensity and the dielectric contrast at optical frequency, 
i.e., 

<p D (r) = -V {I(t)/4\E \% (A6) 

where Vq > equals Vq used for the parametrization of 
the external potential in experiment of Ref. [2^ • 

The strong mutual repulsion of the three constituents 
is mediated via a screened Coulomb interaction 

<pc(r) = -e-* r , (A7) 
r 

where k is the inverse Debye length, r is the center to cen- 
ter distance, and q 2 = [(Z*e ) 2 /(47re er)] e KR ° /(1 + kR s ) 
with Z*eo being the charge on the surface of a single 
particle, e r the dielectric contrast of water, and R s the 
radius of a colloidal particle [43L |44| . Introducing dimen- 
sionless parameters p = r/a and k = na the screened 
Coulomb repulsion scales as <pc(p) = £0 ex P ( — kp)/p, 
where = q 2 /a. 

The internal potential energy of a single trimer, Ec = 
J2i<j Vc (Pij ) and Pij the distance of two colloids within 
the trimer, 

E c /£o= — c~ k ^, (A8) 
P 

is competing with the compression of the external dielec- 
tric potential, Ed = J2i fn{pi) and pt the distance of a 
colloid from the center of mass of the trimer, 

E D /So = -3w (l + 2cos27rp), (A9) 

where vq = Vq/£q. 

The size of a trimer is determined by minimizing the 
total potential energy, Eq + Eh, see Fig. 1211 with respect 
to p. The minimum defines the binding energy Et, which 
for typical experimental parameters exceeds the thermal 
scale by approximately two orders of magnitude - E), ~ 
10 2 k B T, for k = 570 nm, a = 11.5 pm, £ = 10 5 k B T, 
and Vq = 60— llOfc^T. The separation of scales justifies 
the idea to treat the assembly of three colloidal particles 
corresponding to a lattice site as a rigid composite object. 

In the picture of composite objects the trimer-trimer 
interaction enters as residual interaction expected to be 
much smaller than the binding energy. The trimer-trimer 
interaction lifts the orientational degeneracy and is the 
origin of the ordering phenomena under consideration. 
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FIG. 21: Dependence of the binding energy of a trimer on its 
size. The equilibrium trimer extension is determined by the 
minimum of the binding energy, p < 1/3. 



Since the Coulomb interaction is strongly screened it 
is sufficient to consider only neighboring lattice sites. 
For the three different geometrical arrangements of the 
trimers, see Fig. |5J the interaction energy is determined 
by the sum of the nine direct colloid-colloid contributions 
of the rigid triangles, 



E„ 



expf-fcp^s,*]). 



(A10) 



Explicitly, elementary geometrical considerations yield 
for the tables of the distances in dimensionless units 
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where po is the equilibrium size of the trimer. For the 



experimental parameters used above, one obtains E± ~ 
0.01 k B T,E 2 ~ 0.1 k B T,E 3 ~ 1 k B T and, thus, J > 
[see Eq. Q], which shows that phase transitions are to 
be expected in the realistic experimental setups. 

For the case of dimers the chain of arguments is along 
the same lines as above, and one finds for the dimer- 
dimer interactions in the four different geometrical ar- 
rangements exhibited in Fig. 0(b), 
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E » = £oJ2 — r -sexp{-kp li [8,t]), (A12) 

s,t=l P^ S ' l i 
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For the experimental parameters, as used in the exper- 
iment with trimers, one obtains E\ ~ 10 _1 k B T,E 2 ~ 
1 k B T,E 3 ~ 10 k B T,Ei ~ 10 2 fe s T, or equivalently, 
if > and M <C which sets the colloidal dimer system 
in the regime of the herring bone ground state. In order 
to explore the other regions of the phase diagram one 
would have to change the colloid-colloid interaction, for 
example, paramagnetic particles could be used. 

Small residual thermal fluctuations of the compos- 
ite objects, e.g., stretching modes or small oscillations 
around the minimum of a gross orientational state, can 
be absorbed in temperature dependent nearest-neighbor 
interaction energies. For example, one can use a pertur- 
bative approach, as has been done for the elastic con- 
stants in the case of colloidal suspensions in one dimen- 
sional troughs 01 • 
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